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Abstract. We develop an information theoretic interpretation of the number- 
phase complementarity in atomic systems, where phase is treated as a continuous 

■ positive operator valued measure (POVM). The relevant uncertainty principle is 
Oh' obtained as an upper bound on a sum of knowledge of these two observables 

for the case of two-level systems. A tighter bound characterizing the uncertainty 
£^ ' relation is obtained numerically in terms of a weighted knowledge sum involv- 

ing these variables. We point out that complementarity in these systems departs 
from mutual unbiasededness in two signalificant ways: first, the maximum knowl- 
edge of a POVM variable is less than log(dimension) bits; second, surprisingly, 
for higher dimensional systems, the unbiasedness may not be mutual but unidi- 
, rectional in that phase remains unbiased with respect to number states, but not 

^ ■ vice versa. Finally, we study the effect of non-dissipative and dissipative noise on 

' these complementary variables for a single-qubit system. 

t> ' 
oo 

1 Introduction 

Two observables A and B of a ci- level system are called complementary if knowledge of the 
f*"** \ measured value of A implies maximal uncertainty of the measured value of B, and vice versa 

[112] . Complementarity is an aspect of the Heisenberg uncertainty principle, which says that 
for any state ip, the probability distributions obtained by measuring A and B cannot both 
be arbitrarily peaked if A and B are sufficiently non-commuting. Heisenberg uncertainty is 

■ traditionally expressed by the relation 

a- i 

A v ,AA V; i?> -\([A,B})^\, (1) 

where (A^A) 2 = (A 2 )^ — ((A)^) 2 . However, this representation of the Heisenberg uncertainty 
relation has the disadvantage that the right hand side of Eq. (JTJ) is not a fixed lower bound but 
is state dependent. For example, if ip is an eigenstate of A, then both A^A and the right hand 
side of Eq. ([I]) vanish, so that no restriction is imposed on the uncertainty in B. To improve this 
situation, an information theoretic (or "entropic" ) version of the Heisenberg uncertainty rela- 
tionship has been proposed |l|2|3j , which relies on Shannon entropy of measurement outcomes 
as a measure of uncertainty [4I5| . An application of this idea to obtain an entropic uncertainty 
relation for oscillator systems in the Pegg-Barnett scheme [6] has been made in Ref. [7], and 
for entropic uncertainty relations among more than two complementary variables, in Ref. [S]. 

Given two observables A = ^2 a a\a)(a\ and B = ^2f,b\b){b\, let the entropy generated by 
measuring A or B on a state be given by, respectively, H(A) and H(B). The information 
theoretic representation of the Heisenberg uncertainty principle states that H(A) + H(B) > 

2 log ( f(A,B) ) ' wnere f{A, B) = max a ,h | (a\b) |, and H (•) is the Shannon binary entropy. We note 
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that f(A, B) > d 1 ^ 2 , where d is the (finite) dimension of the system. A pair of observables, A 
and B, for which f(A, B) = d~ x l 2 are said to form mutually unbiased bases (MUB) [9110] . Thus, 
any \a) is an equal amplitude superposition in the basis {\b)} and vice versa. Conventionally, two 
Hermitian observables are called complementary only if they are mutually unbiased. Given a 
mutually unbiased pair of Hermitian observables, A and B, the Heisenberg uncertainty relation 
takes the form 

H(A) + H(B)> log d. (2) 

A further advantage of the entropic version of the uncertainty principle over ([1]) is that unlike the 
latter, it is insensitive to eigenvalue relabeling, and depends only on the probability distribution 
obtained by measuring A or B on a given state [3]. 

Even the information theoretic representation ([2]) may not in general be suitable if A or B 
is not discrete, because the continuous analog of H(A), which is H c (p) = — J x dx p(x) log\p(x)], 
is not positive definite, as can be seen from the case where the probability distribution is given 
by p(x) = 2 for x G [0, |] and p(x) — for x G (i, 1], where we find H c (p) = — log 2. It 
is well possible that this pathological behaviour does not afflict classes of physical states of 
interest. In particular, we verified this in the case of the phase distribution of two- and four- 
level atomic systems. However, we are not aware that this is generically true. In any case, this 
potential problem can be generally overcome if the uncertainty principle is expressed in terms of 
relative entropy (also called Kullback-Leiblcr divergence, which is always positive) [11], instead 
of Shannon entropy. An example of where this finds application would be when one of the 
observables, say A, is bounded, and its conjugate B is described not as a Hermitian operator 
but as a continuous-valued POVM. An instance of this kind, considered below in detail, is the 
number and phase of an atomic system. Here we show that the relative entropic definition can 
be used to express complementarity of number and phase, where the notion of complementarity 
is extended to accomodate POVMs. We thus makes this intuitive notion more concrete. Here 
the 'number' variable is analogous to energy in oscillator systems (in the sense of having discrete 
eigenvalues with fixed difference between consecutive values) and amplitude of light field (eg., 
a laser, in the sense of being conjugate to a phase variable). We note that recourse to relative 
entropy is not necessary for a POVM of discrete variables [12], since Shannon entropy is well 
defined in this case. 

The quantum description of phases [13] has a long history |6|14|15|16|17|18] . Pegg and Bar- 
nett [6] , following Dirac |14| , carried out a polar decomposition of the annihilation operator and 
defined a Hermitian phase operator in a finite-dimensional Hilbert space. In their scheme, the 
expectation value of a function of the phase operator is first carried out in a finite-dimensional 
Hilbert space, and then the dimension is taken to the limit of infinity. However, it is not possible 
to interpret this expectation value as that of a function of a Hermitian phase operator in an 
infinite-dimensional Hilbert space |19|20| . To circumvent this problem, the concept of phase 
distribution for the quantum phase has been introduced [19121] , In this scheme, one associates 
a phase distribution to a given state such that the average of a function of the phase operator 
in the state, computed with the phase distribution, reproduces the results of Pegg and Barnctt. 

An interesting question to ask is how mutually unbiased observables behave in the presence 
of noise. Intuitively, one would expect that the uncertainty or entropy of each observable should 
be non-decreasing under the effect of noise. However, this is not generally true, as seen for 
example in the case of a quantum deleter |22I23| , where uncertainty in the computational basis 
vanishes asymptotically during a qubit's dissipative interaction with a vacuum bath. Here we 
study number and phase of atomic systems subjected to both non-dissipative and dissipative 
noise. Noise can be thought of as a manifestation of an open system effect [23]. The total 
Hamiltonian is H = H$ + Hr + H$r , where S stands for the system, R for the reservoir 
and SR for the system-reservoir interaction. The evolution of the system of interest S (in this 
case the atomic system) is studied taking into account the effect of its environment _R, through 
the SR interaction term, making the resulting dynamics non-unitary. The open system effects 
can be broadly classified into non-dissipative, corresponding to the case where [Hs,Hsr] = 
resulting in decoherence without dissipation or dissipative, corresponding to the case where 
[Hs,Hsr] resulting in decoherence along with dissipation [25] , 
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A class of obscrvablcs that may be measured repeatedly with arbitrary precision, with the 
influence of the measurement apparatus on the system being confined strictly to the conju- 
gate observables, is called quantum non-demolition (QND) or back-action evasive observables 
|26|27|28(29] . Such a measurement scheme was originally introduced in the context of the de- 
tection of gravitational waves [30131] . The non-dissipative open system effect described above 
would be a QND effect. Since they describe dephasing without dissipation, a study of phase 
diffusion in such a situation would be important from the context of a number of experimental 
situations. A study of the quantum phase diffusion in a number of QND systems was carried 
out in Ref. [32] using the phase distribution approach. In Ref. [23), the above study was ex- 
tended to include the effect of dissipation on phase diffusion. This would be under the rubric 
of a dissipative open system effect, described above. 

In this paper we study three broad, related problems: first, we formulate a novel character- 
ization of the Hciscnbcrg uncertainty relationship in terms of Kullback-Lcibler divergence (or 
relative entropy). Second, we motivate it by applying it to a study of complementarity in an 
angular momentum system, which involves a continuous variable POVM; lastly, we study the 
behavior of complementary variables when subjected to dissipative and non-dissipative (purely 
dephasing) noise. 

The plan of the paper is as follows. In Section [21 we introduce the concept of phase distri- 
bution in an atomic system which would be used subsequently. In Section [31 we motivate and 
develop an information theoretic representation of complementarity as applied to a two-level 
atomic system, with a brief discussion of a four-level atomic system. Since any system of in- 
terest would, inevitably, be surrounded by an environment which would effect its dynamics, it 
is of relevance to discuss the above ideas of complementarity in the context of open quantum 
systems. We do this in Section [4] by recapitulating relevant work [23 25 32 33 134] on open quan- 
tum systems, of relevance here. Section |4 . 1 1 deals with the non-dissipative open system effect, 
described by the phase damping channel |4|25|32|33] . and Section [4721 discusses the dissipative 
open system effect, described by the squeezed generalized amplitude damping channel [33|34j . 
The reason for the above terminologies is the connection of the dynamics generated by these 
processes with the noise effects pertinent to quantum information |33j . For completeness, we 
relegate some technical details pertaining to these noisy channels to Appendix A and B, where 
the physical processes underlying these channels are also briefly discussed. In Section[5]we make 
our conclusions and discuss some open questions coming out of our work. 

2 Phase distribution 

It is not possible to interpret the expectation value of a function of the phase operator, in 
the Pegg and Barnett scheme [6] , as the expectation value of a function of a Hermitian phase 
operator in an infinite-dimensional Hilbert space |19|20j . This motivates the introduction of the 
phase distribution for oscillator systems [19121] . Interestingly, the concept of phase distribution 
can also be extended to atomic systems [35], which we study here . The phase distribution 
V(4>), <t> being related to the phase of the dipole moment of the system, is given by 



where V(4>) > and is normalized to unity, i.e., d(f>V(4') = 1- I n the above, j is the 
angular momentum of the atomic system. The quantity <j> is important in the context of atomic 
coherences and the interfcrometry based on such coherences [35]. Here Q(9, <fi) is defined as 



where \9, <f>) are the atomic coherent states [36137] given by an expansion over the Wigner-Dickc 
states [38] , which are the simultaneous eigenstates of the angular momentum operators J 2 and 
J z , as 




(3) 



Q(M) = <Mp*IM>, 



(4) 




(5) 
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It can be shown that the angular momentum operators J^^J-q and Jq (obtained by rotating 
the operators J X ,J V and J z through an angle 9 about an axis h = (sin 0, — cos <f>, 0)), being 
mutually non-commuting, obey an uncertainty relationship of the type (J|)(J^) > j(Jq). 
Atomic coherent states (obtained by rotating the Wigncr-Dicke states via 8 and <j> as above) 
are precisely those states that saturate this bound, from which the name is derived |35j . For two 
level systems, they exhaust all pure states, whereas for larger dimensions, this is no longer true. 
Using Eq. (H|) in Eq. ([3]), with insertions of partitions of unity in terms of the Wigner-Dickc 
states, we can write the phase distribution function as [32j 



V ^ = ^TJ d6sin9 E (0^\j,n)(j,n\p s (t)\j,m)(j, m \e,cb). (6) 

^ n 7 m=—j 

The phase distribution V(<fi), taking into account the environmental effects, have been studied 
in detail for QND as well as dissipative systems in [32123] for physically interesting initial 
conditions of the system S, i.e., (a). Wigncr-Dicke state, (b). atomic coherent state and (c). 
atomic squeezed state. 



3 Information theoretic representation of complementarity 

The relative entropy associated with a discrete distribution f(j) with respect to a distribution 
g(j) defined over the same index set, is given by 

s(/llri = £/Cfliog(^). ( 7 ) 

It can be thought of as a measure of 'distance' of distribution / from distribution g in that 
S(f\\g) > 0, where the equality holds if and only if f(j) = g(j) [4]. Consider random variable 
F with probability distribution /. We will define R(F) as the relative entropy of / with respect 
to the uniform distribution -j, i.e., 

R(F) = R[f(j)]=J2f(mog(df(j)). (8) 

3 

As a measure of distance from a uniform distribution, which has maximal entropy, R{F) can be 
interpreted as a measure of knowledge, as against uncertainty, of the random variable described 
by distribution /. The following theorem re-casts Hcisenbcrg uncertainty principle in terms of 
relative entropy. 

Theorem 1 Given two mutually unbiased Hermitian observables A and B, the uncertainty 
relation H|j is equivalent to 

R(A) + R(B) < logd, (9) 
where d is the (finite) dimension of the system. 

Proof. Let the distribution obtained by measuring A and B on a given state be, respectively, 
{pj} and {qk}- The l.h.s is given by 

j k 

= ~[H(A) + H(B)} +21ogd (10) 

^- 2l0g (7(l^)) +2l0gd (n) 
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This is the general result for any two non-commuting Hermitian obscrvablcs. If A and B are 
mutually unbiased, then f(A, B) = d~? , and the theorem follows. It follows from the concavity 
of H, and thus from the convexity of R, that the inequality Eq. © derived for pure states holds 
also for mixed states. ■ 

Physically. Eq. ((9|) expresses the fact that simultaneous knowledge of A and B is bounded 
above by log d. This is in contrast to inequality ((2|) , which is bounded below, being a statement 
on the sum of ignorances or uncertainties. Both are equivalent ways of expressing the fact that 
the probability distributions obtained by measuring A and B on several identical copies of a 
given state cannot both peak simultaneously. 

In terms of R, two Hermitian observables A and B of a d-level system are called mutually 
unbiased if the maximal knowledge of the measured value of A, given by log<i bits, implies 
minimal knowledge of the measured value of B, given by bits, and vice versa. In anticipation 
of the introduction of POVMs instead of Hermitian observables, we will find it convenient to 
weaken the definition of mutual unbiasedness and call two variables A and B (of which one or 
both of them may be a POVM) as quasi-mutually unbiased if the maximal knowledge of the 
measured value of A implies minimal knowledge of the measured value of B, and vice versa. 
The maximum knowledge no longer being logd bits, but lesser, the pair A and B may be called 
quasi-mutually unbiased bases (quasi-MUB's), an extension of the concept of MUB from the 
case of orthonormal bases to that of non-orthonormal bases. 

If two observables are not mutually unbiased, then logd does not bound from above the 
knowledge sum Rt = R{A) + R(B), and there exist states such that the corresponding sum 
satisfies Rt > log d. Intuitively, this is because in the case of two observables that are not mutu- 
ally unbiased, knowledge of the two observables pertaining to a given state may simultaneously 
peak. For example, consider the qubit observables <r z and n ■ a in the Hilbert space C 2 , where 
n = (sin 9, 0, cos 9) and a is the vector of Pauli matrices. It can be seen using Eq. (fTU|) that 
any eigenstate of n • a corresponds to the knowledge sum Rt = 2 — H (cos 2 (9 / '2)) . This sum is 
greater than one, except for 9 = tt/2, which corresponds to the mutually unbiased observable 

Eq. ([7]) has a natural extension to the continuous case, given by 

S(f\\g) = J 4p/(p) log (fflV (12) 

As in the discrete case, wc define R(f) as relative entropy setting g(p) to a continuous constant 
function. In particular, the relative entropy of V((f>) with respect to a uniform distribution ^- 
[32I23| over is given by the functional 

R[V{4>)]= #P(^)log[27rP(0)], (13) 

where the log(-) refers to the binary base. 

We define minimum entropy states with respect to an observable as states that yield the 
minimum entropy when the observable is measured on them. In the context of relative entropy, 
these states can be generalized to what may be called maximum knowledge (MXK) states, which 
are applicable even when the measured variable is continuous. For projector valued measure- 
ments (PVMs), clearly any eigenstate is a MXK state, with a corresponding entropy of zero and 
knowledge R = log<i. PVMs, projectors to the eigenstates of a Hermitian operator represent- 
ing an observable, satisfy three axiomatic requirements: they are positive operators that form a 
partition of unity; further they satisfy the orthonormalcy condition / / '/. = SjkPj, where Pj is a 
measurement operator. The last property implies the idempotency of projectors, which captures 
the idea that projective measurements are repeatablc. From a quantum information perspec- 
tive, it is useful to consider generalized measurements in which the operator elements M m may 
not be orthonormal, but satisfy the completeness condition J^m M^Mm = I and M^M m > 
[4]. In the context of a qubit, for a generalized measurement, the knowledge corresponding to 
a MXK state can be less than 1, i.e., i?(|MXK)) < 1. For a PVM, we have J?(|MXK)) = 1, 
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whereas a POVM considered here is a measurement strategy such that i?(|MXK)) < 1. The 
reason is that whereas PVM is an orthonormal resolution of unity, a POVM forms a non- 
orthonormal resolution of unity |39) . POVMs are useful elsewhere, in quantum information, as 
general measurement strategies for optimally distinguishing states [4] . 

A plot of Rcf, = R[P ((/))] for a two-level atomic system in an atomic coherent state \a',/3') 
with V(<f>) = [l + f sin(a') cos(/3' - 0)] |32l23j . is given by the dashed curve in Figure (fT|). 
We note that R<j, has no dependence on f3' because /3' occurs in V{4>) only as the translation 
<p — f3' , and R^ is translation invariant, i.e., unchanged under the transformation <fi — ► <f> + A. 
The maximum knowledge R^ of about 0.245 occurs at a 1 = n/2. The corresponding continuous 
family of states |7r/2,/?') forms the MXK states or quasi- eigenstates of the phase observable. 
These are equatorial states on the Bloch sphere, having the form -^(|0) + e 4 *°|l)). That R$ is 
less that 1 for these states reflects the fact that here phase <j> is a POVM. 

In analogy with the oscillator case, the Wigner-Dicke or excitation states may be thought 
of as 'number states', thereby making J z the 'number observable', whose distribution is p(m), 
given in Eq. (|14[) . The 'number' distribution given by 

p(m) = (j,m\p s (t)\j,m), (14) 

is considered as complementary to T > {4>) [35 . It is of interest to ask whether they are comple- 
mentary in the sense of MUBs. 

In the manner of Eq. |(5J), we can define R m = R[p(m)] as knowledge of the number variable. 
We note that J z and phase <j> have a reciprocal behavior reminiscent of MUBs: the eigenstates 
of J z , i.e., Wigner-Dicke states, correspond to minimal knowledge R<p(= 0), as seen from the 
dashed curve in Figure ([1]). This can be seen by noting that for the Wigner-Dicke states \j, rh), 
the phase distribution is [32] 

m = ^ ( ? l^) B [j + m+ l,i - m + 1] = ± (15) 

where B stands for the Beta function. Thus, it follows via Eq. (|13[) that the knowledge R<j, 
vanishes. Conversely, we note that the states which minimize R<j> are the Wigner-Dicke states. 
To see this, we observe that if V((f>) is constant, then in Eq. ©, each term in the summation, 
which is proportional to e , % \ m — n )4>^ mu st individually be independent of <p. Since <f> is arbitrary, 
this is possible only if m = n, i.e., the state p s is diagonal in the Wigner-Dicke basis. Thus, 
MXK states of m correspond precisely to minimum knowledge (MNK) states of <j>. 

The plot of relative entropy R m for all atomic coherent states is given by the bold curve in 
Figure |T]). The equatorial states on the Bloch sphere, the MXK of 0, are precisely equivalent to 
the MNK states of m (characterized by R m = 0), as can be seen from comparing the dashed and 
bold curves in Figure |T]). Thus number and phase share with MUBs the reciprocal property 
that maximum knowledge of one of them is simultaneous with minimal knowledge of the other, 
but differs from MUBs in that the maximum possible knowledge of <f> is less than log(d) = 1 
bit, essentially on account of its POVM nature. 

Two variables form a quasi-MUB if any MXK state of either variable is an MNK state of 
the other, where the knowledge of the MXK state may be less than \ogd bits. Thus, J z and 4> 
are quasi-MUB's (but not MUB's), and are complementary in the extended sense. 

From the dot-dashed curve in this Figure, we numerically find an expression of the uncer- 
tainty principle to be 

R T = R<p + R m < 1 (16) 

for all states (pure or in general mixed) in C 2 , in analogy with Eq. ([9]). The inequality is 
saturated for the Wigner-Dicke states. 

As an expression of the uncertainty principle, the relation (|16p still leaves some room for 
improvement. First, it is not a tight bound. In particular, for equatorial states it permits i?^ 
to be as high as 1 , whereas as seen from the dashed curve in Figure [lj the maximum value 
of Rtf, f» 0.245. We note that the bound cannot be tightened simply by decreasing the r.h.s, 
since it is saturated for Wigner-Dicke states. Further, the variable takes values in the interval 
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Fig. 1. Plot of knowledge R with respect to a' of a qubit starting in an atomic coherent state \a',j3'). 
Knowledge R is symmetric with respect to (3' , which is therefore not depicted in the figure. The 
individual curves are: R^, = R[P(<f>)] (dashed curve), R m = R[P(m)] (bold curve), the knowledge sum 
R4, + Rm (dot-dashed) and the weighted knowledge sum Rs(n) = fiRcf, + Rm with [i = 4.085 (dotted 
curve). The maximum of R$ is not 1 but 0.245 bits, which occurs at a' = n/2, i.e., the equatorial states 
on the Bloch sphere, which are thus the MXK states or quasi-eigenstates of cj>. The minimum value of 
Rtf, = occurs at a' = and tt, corresponding to the Wigner-Dicke states. 



[0, 2tt] irrespective of the dimensionality of the Hilbcrt space, unlike m, which takes d values. 
Consequently, , unlike R m , is not seen to be bounded by the dimension of the Hilbcrt space 
in a straightforward way. To see that in general R[p(x)] increases without bound, consider the 
probability distribution p(x) = xo > 1 in x G [0, — 1 and p(x) = in x G ( — , 1], for which we 
find R(p(x)\\\) = logxo- 

One way to address these problems is to generalize (|16p to a family of inequalities, parametrized 
by /i > 0, of the form 

R s (ji) = ^ + Rm<l (17) 

for all states in C 2 . We find that the largest value of fi such that inequality (|17[) is satified 
over all state space is (r < f > )~ 1 « 4.085, where is the value of R^ for the equatorial states, the 
MXK states of 4>. A plot of Rs{l/r^) over pure states is shown as the dotted curve in Figure 
[TJ Comparing this with the dot-dashed curve in Figure [TJ we find that Rs(fJ-) is bounded more 
tightly than R T = R s (l). 

From Figure [TJ we find that the two Wigner-Dicke states and all equatorial states may 
be regarded as coherent with respect to the number-phase pair, in that they maximize the 
knowledge sum and are thus closest to classical states. We note of course that this definition 
of state coherence differs from the conventional one for atomic states, defined with respect to 
angular momentum operators. Unless we use p-R^ in place of R^, only the Wigner-Dicke states 
could be called coherent in the new sense. 

We now briefly extend the entropic version of complementarity to a higher spin system, 
which is seen to present a new feature. We consider a spin-3/2 (four-level) system, whose 
general state is given by the ansatz 

W=rae*-||,-§>+r^|| 1 -i>+r T e^|| > +i)+r,||,+|>) (18) 
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where r\ + r \ + rl + rf = 1, and a global phase is omitted. Using Eq. (fT8|) in Eq. ((6]), we find 



1 


'1 


/ 157rr Q r,3\ 


7T 


2 + 


V 32V3 / 



C os(0 - e a 6p) + y-^=-j c 08 ^ - 



^ cos(30 - Q ) + [-^j cos(0 - ^ + 7 ) + I ) cos(20 - Bp) 



32\/3 



'7 



(19) 



As before, we compute 'number' knowledge i? m (r Q , rg, r 7 ) using Eq. ([8]), and phase knowl- 
edge R^(r a ^r/3,rj^9 a ,9p^9y) using Eq. (|13[) . It may be verified that for 'number' states (for 
which r a or 773 or r 7 or r,5 is 1), R^, = 0. In fact, it may be seen from Eqs. |5"]). ([S]) and 
(fTS"]) . that a general property of atomic systems is that a Wigner-Dicke state is equivalent to 
a MNK phase state in any finite dimension. On the other hand, numerically searching over all 
possible states of the form (fTS)) . we find that the maximum value 0.86 bits of R<f, occurs at 
l?(r Q = 0.36,773 = 0.61, r 7 = 0.61, Q = tt, 6p = 0, 6* 7 = n), which is not an equal amplitude 
superposition of 'number' states. Thus, remarkably, for the spin-3/2 case, MXK phase states 
do not correspond to MNK 'number' states, even though the converse is true. We expect that 
this unidirectionally (as against mutually) unbiased behavior will persist even for higher spin 
systems. Phase and 'number' therefore do not here form a quasi-MUB as defined for the single 
qubit case, and may be considered complementary only in an even more weak sense. This is in 
contrast to the case where the observables are Hermitian, where five MUBs are known to exist 
in four dimensions |10| . 

As in the two-level case, one way to address this problem is to generalize ([To]) to a family 
of inequalities, parametrized by fj.2 > 0, of the form 

Rs(»2) = lizR* + Rm < 2 (20) 

over all states in C 4 . Our strategy is to numerically search over all states of the form (fT5)) - other 
than the Wiger-Dicke states, where R m = 2 and R^ = 0- in order to determine the largest 
value of fi2 such that inequality ([20|) is just satificd, i.e., the inequality must be satisfied at all 
points, with the equality being valid for at least one point. By this method, we find /12 = 1.973 
with the maximum Rs(n2) of 2 occuring at tp p = ^{r a = 0.24, rp = 0.64, r 1 = 0.68, 9 a = 
n, 9p = 0, 6* 7 = n). As states that maximize the knowledge sum Rs(fi2), we may regard ip p and 
the Wigner-Dicke states as coherent states from the viewpoint of number-phase entropy. 



4 Application to open systems 



Here we study the effect of noise coming from open quantum system effects, on the atomic 
number-phase complementarity developed in the previous section. The noise effects we consider 
come from non-dissipative as well as dissipative interactions of the atomic system S with its 
environment which is modelled as a bath of harmonic oscillators starting in a squeezed thermal 
state [25 33 34 . This enables us to consider the effect of bath squeezing on the complementarity. 
We briefly recapitulate previous work [23 25 32 33 3j[] related to the effect of various noisy 
channels on the 'number' and phase distributions. In Section |4~T1 we consider the effect of the 
phase damping channel which is the information theoretic analogue of the non-dissipative open 
system effect [25133] while in Section 14.21 we consider the effect of the squeezed generalized 
amplitude damping channel which is the information theoretic analogue of the dissipative open 
system effect |33|34j . Intuitively, one would expect that open system effects, like measurements, 
cannot increase the knowledge sum. Interestingly, we find that this is not true for certain regimes 
of the squeezed generalized amplitude damping channel. 
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4.1 Phase damping channel 

The 'number' and phase distributions for a qubit starting from an atomic coherent state |a', /?'), 
and subjected to a phase damping channel due to its interaction with a squeezed thermal bath, 
are [28125182] 

p(m) = L^Jj (sin(a72)) 2 W+™)(cos(a'/2)) 2 ^- m ) 



1 



2tt 



-p(<p) = —l+l sin(a') cos(/3' + tot - cf>)e~^^ 



4 



(21) 



R<f, (Eq. I|13p) is invariant under the translation of — ► <f> + a. Setting a = —j3' — ut, we 
find that R<f, is independent of /?'. A derivation of Eq. (j2Tj) can be found in Refs. [32l23j . For 
completeness, the expression for j(t) in Eq. (|2ip is given in Appendix [Al and the physical process 
underlying the phase damping channel discussed. 

Figure [2] depicts the effect of phase damping noise on the knowledge sum Rs- Comparing it 
with the noiseless case (dotted curve in Figure [l}, we find a reduction in the total knowledge 
Rs, as expected. It follows from Eq. (f2"Tj) that R m remains unaffected under the action of this 
channel. Thus, the effect of noise on Rs is due entirely to its effect on i?^, which decreases in 
the presence of noise for all pure states (because j3' does not play any role and because the plot 
represents all possible values of a'). 

Figure [2] shows that squeezing has the detrimental effect of impairing phase knowledge for 
all regimes of the parameter space. This is in marked contrast to the case of the squeezed 
generalized amplitude damping noise, discussed in Section (|4.2[) . Thus, squeezing, like temper- 
ature, has the overall detrimental effect of impairing Rs- This is consistent for the case of a 
QND interaction (which generates a phase damping channel |25|33j ) of the system with its 
environment, i.e., [Hs, Hsr] = 0, as also corrorborated by the observation that squeezing and 
temperature concurrently impair geometric phase [33] and phase diffusion [32|23j and suggests 
that squeezing, like temperature, should adversely affect channel capacity for phase damping 
noise. 



4.2 Squeezed generalized amplitude damping channel 



The 'number' and phase distributions for a qubit starting from an atomic coherent state \a', (3'}, 
and subjected to a squeezed generalized amplitude damping channel [34] due to its interaction 
with a squeezed thermal bath, are [2"3"] . 



p(m= 1/2, i) 
and 

v{4>)- 



i 

2^ 



1 - 



To 



1 



sinV/2) + ^ (l - e" 7 "*) cosV/2), 



7T f 

+ — sin(a') jacosh(ai) cos(0 — j3 ) + wsinh(erf) sin(0 — (3 ) 



-7oxsinh(o;t) cos(<P + + 0) je 



(22) 



(23) 



A derivation of Eqs. (|2"2"|) and (|2"3"1) can be found in Refs. [53]. For completeness, the parameters 
appearing in these equations are given in Appendix [B] where a brief discussion of the physical 
process behind the squeezed generalized amplitude damping channel is also made. 

Figures [3Ja) and (b) depict the effect of squeezed generalized amplitude damping noise 
on iiRrf, (Eq. (j!7|l ). without and with bath squeezing, respectively. Comparing them with the 
noiseless case of Figure [1] (which, it may be noted, is unsealed by /i), we find as expected 
that noise impairs phase knowledge. However, comparing Figure [Hb) with (a), we find that 
squeezing has the beneficial effect of relatively improving phase knowledge for certain regimes 
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Fig. 2. Plot of Rs = fiRg) + Rm for a qubit starting from an atomic coherent state \a',/3' = 7r/4), 
subjected to a phase damping channel, with temperature (in units where h = ItB = 1) T = 2, 70 = 0.025, 
u>c = 100 and lj — 1.0, with respect to bath exposure time and bath squeezing parameters (Eq. (13611 ) 
r = 1 and a = 0.0. For the phase damping channel, R m remains invariant. Both R m and R$, and hence 
Rs, are independent of 



of the parameter space, and the detrimental effect of relatively impairing them in others. This 
property can be shown to improve the classical channel capacity [34]. Further, bath squeezing 
is seen to render dependent on f3', because, as evident from Eq. (|23|) . (3' no longer appears 
as a translation in (f> when the squeezing parameter \ (Eq. (|39p ) is non- vanishing. On the other 
hand, it follows from Eq. (|22[) that R m is independent of /?', so that Rs is dependent on /?'. 
This stands in contrast to that of the phase damping channel, where inspite of squeezing, Rs 
remains independent of f3' and, furthermore, squeezing impairs knowledge of <j) m regimes 
of the parameter space. 

A point worth noting is that, in contrast to the phase damping channel, in a squeezed gen- 
eralized amplitude damping channel, R m and Rs are not necessarily non-increasing functions 
of time. Figure HJa) depicts the effect of squeezed generalized amplitude damping channel on 
Rs, by bringing out the behavior of Rs as a function of bath exposure time. The dashed curve 
shows that squeezing has a detrimental effect on the knowledge sum, as one would usually 
expect. A surprising departure from this behavior may be noted for the case of the bold curve, 
which corresponds to the action of a dissipativc interaction with an unsqueezed vacuum bath, 
where the knowledge sum Rs increases to 1. This counterintuitive behavior is due to the quan- 
tum deleting action, a contractive map whereby any initial state, including a mixed state, is 
asymptotically prepared in the pure state |i, — i) for vanishing temperature, and a mixture of 
|0) and |1) states for finite temperature, where the asymptotic mixture is determined purely by 
the environmental parameters of T and r, and not by the system's initial state |22| . A similar 
effect was noted in [40j . where in a study of quantum state diffusion of an open system it was 
shown that for a specific noise, due to a particular system-reservoir interaction, there can be a 
reduction in the quantum dispersion entropy leading to localization. 

It follows from the complementarity relation Eq. dTTJ) , that in the asymptotic limit of the 
deleting action, R^ goes to for both |0) and |1), and hence also, by the convexity of R, for 
any mixture that is diagonal in this basis. More generally, it is seen from Figure [UJb) that for 
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Fig. 3. Plot of jiRj, (scaled R4, [Eq. (|17[l]) for a qubit starting from an atomic coherent state \a ,/3 ), 
and subjected to a squeezed generalized amplitude damping channel with temperature (in units where 
h = ks = 1) T — 300, 70 = 0.01, bath exposure time t = 0.1, and lo = 1.0. Figure (a) is for the case of 
zero bath squeezing, and (b) for the case of bath squeezing parameters (Eq. (|39j) r = 1 and $ = it/8. 
Squeezing has the effect of breaking translation symmetry in (3' and improving phase knowledge (i.e., 
reducing phase uncertainty) for certain regimes of the parameter space. 



0.2 




Fig. 4. (a) Plot of Rs = ^R^, + R m with respect to time for a qubit starting from an atomic coherent 
state \oc — n/4,/3' = 7r/4) and subjected to a squeezed generalized amplitude damping channel. The 
temperature (in units where h = ks = 1) T = 0, 70 = 0.025, and U) — 1.0. The bath squeezing 
parameters (Eq. ([39j ) are ^ = 0, and r = (r = 0.5) for the case of bold (dashed) curves, (b) Plot of 
/ii?^ for all a' as a function of time for a qubit starting from a coherent state |a',7r/2) and subjected 
a squeezed amplitude damping channel with the same parameters as above, with T — <P = 0, r = 0.5 
and 7 = 0.05. The i?^,-decreasing effect of increasing T or r is qualitatively the same. 



all initial pure states, R<p falls monotonically. This is to be expected since this noise prepares 
an asymptotic state that lies on the z-axis of the Bloch sphere, which implies by the convexity 
property of R$ and the fact that = for the north and south pole states, that asymptotically 
R<j> = for all initial pure states. 
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5 Conclusions 

In this work, we have investigated the number-phase complementarity in atomic systems from 
an entropic perspective through the number and phase distributions. Here number distribution 
refers to the probability distribution of measurement in the Wigner-Dicke basis (Eq. (|14j)). 
while phase distribution is defined by Eq. We derive an uncertainty principle in terms of 
the Kullback-Leibler or relative entropy R of number and phase with respect to a uniform 
distribution. Since R can be regarded as a measure of knowledge of a random variable, the 
entropic uncertainty principle takes the form of an upper bound on the sum of number knowl- 
edge (R m ) and phase knowledge (i?</>). The choice of relative entropy over Shannon entropy 
was motivated by the fact that the latter is not strictly positive when applied to continuous 
probability distributions. 

In the single-qubit case, number and phase are regarded as quasi-MUBs in the sense that 
any state maximizing knowledge of one variable simultaneously minimizes knowledge of the 
other, but maximum phase knowledge is strictly less than 1 bit (and less than logd bits in d 
dimensions). 

Since R^ is strictly less than one, the relative entropic formulation of the uncertainty prin- 
ciple does not tightly constrain R m . We define a family of inequalities, parametrized by [i (Eq. 
OH)), that improves the upper bound on R m . When \i = 1, we obtain Eq. (|16[) . and get the 
tightest bound for equatorial states (with the right hand side saturated) when /i w 4.085. We 
briefly study the extension of the above concepts to a four-level system, where we find that the 
sense in which number and phase are said to be complementary must be further weakened to 
include unidirectional (but not mutual) unbiascdness. In particular, whereas phase is unbiased 
with respect to number, the converse is not true. 

Finally, we study the complementary behavior of number and phase of a qubit subjected 
to the influence of its environment. For a qubit starting from an atomic coherent state |a',/3'), 
the translation symmetry of R^ in /3 is broken by the introduction of squeezing in the bath, 
for the case of a dissipative system-bath interaction (Figure [31(b) ) , but not in the case of a 
non-dissipativc interaction. In the case of a purely decohering interaction, characterized by a 
phase damping channel, we find that noise invariably impairs the knowledge sum for these 
complementary variables (Eq. (JTTJ) ) , as expected. However, in the case of a squeezed general- 
ized amplitude channel, the knowledge sum can increase in certain regimes. As a particularly 
dramatic illustration, when an initially maximally mixed state \{\\, h\ + ||i — l)(l> — h\) 
is subjected to an unsqueezed vacuum bath, Rs rises from to 1 asymptotically. 

These results could be potentially useful for applications in quantum communication and 
quantum cryptography |41j involving atomic systems. The present work brings forth a number of 
open questions concerning an information theoretic study of complementarity in atomic systems 
involving continuous-valued POVMs, of which we list some here. Of immediate interest is the 
question whether the Shannon entropy of V(4>) remains positive for all possible pure and mixed 
states. If yes, then one may revert back from the use of the knowledge variable R to that of 
entropy S. Also of interest is to analytically derive the bounds on the weighted knowledge sum, 
which we have obtained here numerically. Finally, it is of interest to explore the full scope and 
implications of one-way biasedness, and its connection to complementarity. 



A Phase damping channel 

Consider the Hamiltonian 



Here Hs, Hr and Hsr stand for the Hamiltonians of the system, reservoir and system-reservoir 
interaction, respectively. Hs is a generic system Hamiltonian which can be specified depending 



H = Hs + Hr + Hsr 




(24) 
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on the physical situation. b k , b k denote the creation and annihilation operators for the reservoir 
oscillator of frequency oj k , <7fc stands for the coupling constant (assumed real) for the interaction 
of the oscillator field with the system. The last term on the right-hand side of Eq. (1) is a 
renormalization inducing 'counter term'. Since [Hs,Hsr] = 0, the Hamiltonian (1) is of QND 
type. The system-plus-reservoir composite is closed and hence obeys a unitary evolution given 

by 



P(t) 



-iHt/h 



p(0)e 



iHt/h 



where 



p(0) = P s (0) P r(0), 



(25) 



(26) 



i.e., we assume separable initial conditions. Here Pr(0) is the initial density matrix of the 
reservoir which we take to be a squeezed thermal bath given by 



p R (0)=S(r,*)pt h &(r,$), 



(27) 



where 



Pth 



(28) 



is the density matrix of the thermal bath at temperature T, with (3 = l/(fcsT), ks being the 
Boltzmann constant, and 



S(r k ,$ k ) = exp 



(29) 



is the squeezing operator with r k , <P k being the squeezing parameters |42j . Here we take the 
system to be a two-level atomic system, with the Hamiltonian 



H, 



(30) 



a z being the usual Pauli matrix. The reduced density matrix of the system, in the basis of the 
Wigncr-Dicke states \j, m) , after time t is [33] 



cos 2 (^) 

Pm ' n{ ' 1 I S m(0 o ) e 'M+*o) e -(M 2 7(t) 



1 sm(6 )e~ i( - UJt+ ' 1 "^ e - ^^ 1 ^ 
h 

2 



sin 2 (^) 



(31) 



from which the Bloch vectors can be extracted to yield 



(a x (t)) = sin(0o) cos(u;f + Me^^, 
(cr y (t)) = sin(e ) sin^t + ^K^ 7 ^ 
(a z {t)) =cos(# ). 



(32) 



Here jit) comes due to the interaction with the environment and for the case of an Ohmic bath 
with spectral density 



(33) 



where 70 and uj c are two bath parameters characterizing the quantum noise, it can shown that 
using Eq. (f3"3"|) one can obtain [35] in the T = limit, 



7 (t) = ^-cosh(2r) ln(l+w 2 i 2 )-— sinh(2r) In 

27T 47T 



(l + 4c e 2 (t-a) 2 ) 
(l + ^(t-2a)2) 2 



-— sinh(2r) ln(l+4a 2 w 

47T 



2 -c 2 ), 



(34) 
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where the resulting integrals are defined only for f > 2a. In the high T limit, 7(f) can be shown 
to be [25] 



7(f) = — cosh(2r) 



2oj c t tan 1 (oj c t ) + In 



1 



1 + w 2 f 2 



— sinn (2r J 



4w c (f — o) tan 1 (2u; c (f — a)) — 4w c (f — 2a) tan 1 (w c (f — 2a)) + 4aw c tan 1 (2aw c ) 



In 



[l+cj 2 (f-2a) 2 ' 
[l + 4cj2(f- a )2 



In 



1 



1 + 4a 2 w 2 



(35) 



where, again, the resulting integrals are defined for f > 2a. Here we have for simplicity taken 
the squeezed bath parameters as 



cosh(2r(cj)) = cosh(2r), sinh (2r(u)) = sinh(2r), 
<P(tu) = auj. 



(36) 



where a is a constant depending upon the squeezed bath. The results pertaining to a thermal 
bath can be obtained from the above equations by setting the squeezing parameters r and <P 
to zero. o x , a y , a z arc the standard Pauli matrices. It can be easily seen from the above Bloch 
vector equations that the QND evolution causes a coplanar, fixed by the polar angle 8q, in-spiral 
towards the z-axis of the Bloch sphere. This is the characteristic of a phase-damping channel 

HI- 



B Squeezed generalized amplitude damping channel 



Here the reduced dynamics of the two level atomic system (|30[) interacting with a squeezed 
thermal bath under a weak Born-Markov and rotating wave approximation is studied. This 
implies that here the system interacts with its environment via a non-QND interaction, i.e., 
[Hs,Hsr] 7^ such that along with a loss in phase information, energy dissipation also takes 
place. The evolution has a Lindblad form which in the interaction picture is given by [24] 

j t P s {t) = 7o (JV + l) (a-ff(t)* + - \<J+<J- P s {t) - \p s {t)a + a. 

+ la N (<j + p s {t)<j- - ^a-a + p s (t) - ^p s {t)a-a+ 
- l0 Ma +P s (t)a+ - 7 MV_ / /(f)a_. (37) 



Here 



and 



TV = 7V t h(cosh 2 r + sinh 2 r) + sinh 2 r, (38) 
M = -I sinh(2r)e i *(2JV tfc + 1), (39) 



^= e>to/(fcflr) _ 1 , (40) 



where Nth is the Planck distribution giving the number of thermal photons at the frequency w, 
and r, are squeezing parameters of the bath. The case of a thermal bath without squeezing 
can be obtained from the above expressions by setting these squeezing parameters to zero. 70 
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is a constant typically denoting the system-environment coupling strength. This equation can 
be expressed in a manifestly Lindblad form as 

, 2 

j/tt) = E ( 2R iP SR ] ~ R \ R jP S ~ P SR ) r j) > (41) 
j=i 

where Ri = (7o(iVth + l)/2) 1 / 2 i?, R 2 = (7 A r th/2) 1/2 i? t . Here R = ct_ cosh(r) + e i,p a + sinh(r), 
and a± = \ (a x ± i<J y ). If T = 0, so that _/V t h = 0, then i?2 vanishes, and a single Lindblad 
operator suffices. The fact that the above equation can be expressed in the form (|4T|) guarantees 
a Kraus or operator- sum representation [4] for the evolution of the reduced density matrix. It 
can be seen that the reduced density matrix, obtained by solving Eq. (157)) in the Bloch form, 
shrinks towards the asymptotic equilibrium state p aS ymp, given by 

P =( 1 ~ P °] (42) 

Pasymp — I q I j 1^,1 

where p = | 1 + ( 2 iv+i) ' ^ or casc 01 zero squeezing and zero temperature, this action 
corresponds to an amplitude-damping channel [4133] with the Bloch sphere shrinking to a point 
representing the state |0) (the south pole of the Bloch sphere) while for the case of finite T 
but zero squeezing, the above action corresponds to a generalized amplitude-damping channel 
[4|33j with the Bloch sphere shrinking to a point along the line joining the south pole to the 
center of the Bloch sphere. The center of the Bloch sphere is reached in the limit of infinite 
temperature. Thus, the interaction with the environment provides a contractive map, such that 
the asymptotic state is pure (p = 1), corresponding to the deletion action [22], or mixed (p < 1), 
depending on environmental conditions. For finite T and bath squeezing, the above corresponds 
to a squeezed generalized amplitude damping channel |34j . 
In Eq. (|23p . the parameter a is given by 



a = y^lMp-^, (43) 

while 

7^ = 7 o(2iV+l). (44) 
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